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ABSTRACT: We prove a strong form of the equivalence of ensembles for the invariant 
measures of zero range processes conditioned to a supercritical density of particles. It 
is known that in this case there is a single site that accomodates a macroscopically 
large number of the particles in the system. We show that in the thermodynamic limit 
the rest of the sites have joint distribution equal to the grand canonical measure at 
the critical density. This improves the result of Grofikinsky, Schutz and Spohn, where 
convergence is obtained for the finite dimensional marginals. We obtain as corollaries 
limit theorems for the order statistics of the components and for the fluctuations of the 
bulk. 
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1 Introduction 

In a landmark paper of 1970, Spitzer |Sp| introduced five particle systems undergoing simple in- 
teractions and initiated a research project to rigorously analyse their equilibrium and dynamical 
properties. One of the systems he proposed was the zero range process, a model in which particles 
leave any given site at a rate g(k) that only depends on the number k of particles present at the site, 
hence the name. The attention was initially drawn to the existence of the dynamics under general 
conditions, the identification of invariant measures and the establishment of the hydrodynamic limit. 
All these questions have been successfully addressed, at least in the attractive case when the rate 
function g(-) is increasing. A comprehensive review of these results can be found in [KLj . 

Over the last decade, there has been an increasing interest in zero range processes such that the 
rate g(-) decreases with the number of particles. This can be thought as introducing a mechanism 
of effective attraction between the particles, that if strong enough, i.e. when the rates decrease 
sufficiently fast, can lead to phenomena of condensation - a transition to a phase where a single site 
contains a finite fraction of the particles in the system. This type of condensation appears in diverse 
contexts such as traffic jamming, gelation in networks, or wealth condensation in macroeconomies, 
and zero range processes or simple variants have been used as prototype models. Evans and Hanney 
EH| provide an excellent review on this subject. 



A phase transition in this class of zero range processes can be already observed at the level of 
the invariant states. It is known [JMP1 IGSSj that when the density of particles exceeds a critical 
value p c , the invariant measures of the process concentrate on configurations where a macroscopic 
proportion of the total number of particles forms a randomly located cluster. In this article wc 
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analyse the thermodynamic limit of the invariant measures of the process conditioned to having a 
supercritical density, that is we let the number of sites L and the number of particles N grow to 
infinity in such a way that N/L — > p > p c . 

Given the particle and site numbers N and L as above, the invariant state of the process is identified 
as the product of L copies of a measure v$ c supported on the integers, conditioned to adding up to 
N. When the particle density N/L is higher than p c we are conditioning on an atypical event, and 
the problem can be described as Gibbs conditioning for a measure having no exponential moments. 
GroBkinsky, Schutz and Spohn [GSS identified the typical configuration of a finite subsystem by 
proving an equivalence of ensembles property. Remarkably, the effect of the conditioning on the 
finite subsystem disappears in the thermodynamic limit. This happens because the rare event is 
most likely realised by a large deviation of the maximum component. A similar result was proved by 
Ferrari, Landim and Sisko [FLS when the number of sites is fixed while the particle number grows 
to infinity, and by GroBkinsky [G] for systems with two particle species. 

The fact that convergence to a product measure holds for the finite dimensional marginals is stan- 
dard when the equivalence of ensembles or the Gibbs conditioning principle are satisfied. It is crucial 
that the size of the subsystem amounts to a vanishing fraction of the whole. Indeed, the result often 
fails to hold when this is not the case (cf. Proposition 2.12 in DZ ). The main result in this article 
is an unusually strong form of the equivalence of ensembles. Precisely, we prove (Theorem [T]) that 
in supercritical zero range processes the effect of conditioning is entirely absorbed by the maximum 
component, in the sense that the joint distribution of the remaining sites converges to a product 
measure. We then derive some interesting corollaries from this result. 

This distinctive behavior can be attributed to the fact that the marginals v < j >c of the unconditional 
distribution are subexponential. Indeed, the proof of Theorem [T] relies on a Local Limit Theorem in 
the form of equation ([8]), a result that requires little more than subexponentiality. 

2 Notation and results 

Zero range processes are interacting particle systems evolving on a set of sites A. Particles perform 
random walks on A interacting only with particles sitting on the same site through the following 
rule: the rate at which a particle leaves a site depends on the number of particles at that site. Given 
a function g : No = {0, 1, 2, . . .} i— > R+ and a transition probability p(-, •) on A x A, the dynamics of 
the process can be described as follows. If there are k particles at a site x, then independently of the 
configuration on the other sites, a particle leaves x after an exponential waiting time with rate g{k). 
A target site is chosen according to p(x, •), the particle jumps there and the process starts afresh. 

A zero range process can be rigorously defined as a Markov process on the state space Xa = Nq. A 
point r\ in Xa can be thought of as a configuration of particles on A, with r\ x denoting the number 
of particles at the site i£A. Regarding the jump rate function g(-) and the transition probabilities 
p(-, •), we assume that 

g : N Q h-> R + is such that g(k) = k = 0, 

and 

p : A x A i— > [0, 1] is such that p(x, y) = p(y, x) = 1, Vie A. 

j/GA y€A 

In order to avoid degeneracies we further assume that the random walk on A with transition proba- 
bilities p(-, •) is irreducible. In this article we only consider finite sets A, in which case we can define 
a process starting from any initial configuration rj £ Xa- 
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The infinitesimal generator of the zero range process is then given by 



where 



Zero range processes possess a family of invariant product measures with site marginals given by 

1 <b k 




. [rix = k] = 



Z{<t>)g{k)\ ' 



where g(k)\ = Yl m =i #( m )- Each of these measures is usually referred to as the grand-canonical en- 
semble corresponding to the fugacity (f>, and they can be defined for any c/> in the range of convergence 
of the power series 



The expected number of particles per site is given by 

p(<f>) = E"* [ Vx ] 



^ OO 



It can be easily verified that p is a strictly increasing function of 4>. 

Let 4> c < +00 denote the radius of convergence of Z(<fi). If Z(cf> c ) :— lim^j,^ Z(4>) — oo, it can 
be proved [KLj that p c := lim^|^ t , p((f>) = oo. If on the other hand Z{<f> c ) is finite, it is possible that 
p c is also finite. In this case none of the grand-canonical measures corresponds to a particle density 
higher than the critical p c , and the system undergoes a phase transition |JMP[[(?55] from a fluid to 
a condensed phase, in a sense to be made precise later. 

To fix ideas, we consider here a reference model such that both Z((f> c ) and p c are finite that was 
originally proposed by Evans [E] . In the last section we discuss how our results apply to a number 
of other systems with finite critical density. 

In Evans' model the jump rates are given by 

9( k ) = S n T , n CO 



With this choice of g, one gets 

T(b + k+ 1) 



9(k) 



r(6 + i)fc! r(6 + i) : 



if r(-) denotes the standard Gamma function. The critical fugacity (f> c is equal to 1, the parti- 
tion function Z(4>) is finite at (f> c if b > 1, and the critical density p c is finite if b > 2. Since we 
are interested in systems with finite critical density we will assume throughout this article that b > 2. 
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Due to the conservation of the number of particles by the dynamics, the state space is partitioned 
into finite invariant subspaces, where 



Sl(v) 



is constant: Xa,jv = {v <= : Sl(t}) — N}. On each of these subspaces the zero range process 
is irreducible and has a unique invariant measure which we denote by p N,L . We will refer to the 
measures p N ' L as the canonical ensembles. They can be explicitly computed, but they can also be 
obtained by conditioning the grand-canonical ensembles on the total number of particles. That is 



.N.L 



[•]=*£[• | S L (V)=N]. 



Note that the right hand side of the last equation does not actually depend on <j>. A natural object 
of interest is the behavior of these measures in the thermodynamic limit, as N, L — > oo in such a 
way that the average particle density N/L converges to a constant p. 

When p < p c there exists a fugacity 4> such that p = p((j>) and the standard equivalence of en- 
sembles for independent random variables holds [KLJ . That is, the finite dimensional marginals of 
the canonical ensembles p N ' L converge to the grand-canonical ensemble corresponding to fugacity 
4>. The equivalence of ensembles for (super)critical densities (p > p c ) was established by Grofikin- 
sky, Schiitz and Spohn [GSS . Using relative entropy methods they prove convergence of the finite 
dimensional marginals of p, * to the grand-canonical ensemble at critical fugacity. 



Furthermore, it has been proved [JMPilGSSllG] that when the density is supercritical a condensation 
phenomenon emerges. Precisely, if p > p c and e > then 



lim u 



N.L 



— max t) x > p - p c 

L x£A 



(2) 



This is to be contrasted with the size of the largest component in the case below criticality, which is 
of order log(L) | JMPj . The comparison gives a precise meaning to the phase transition experienced 
by the system, and is reminiscent of the Erdos-Renyi results on the largest cluster of a random graph. 



The heuristic picture suggests that at supercritical densities the bulk of the sites is distributed 
according to independent copies of i^ , while a single randomly located site accumulates all the 
excess mass. The results mentioned above do not fully justify this picture however, because conver- 
gence to the grand-canonical ensembles is only obtained at the level of finite dimensional marginals. 
Hence, questions that require knowledge of the full limiting distribution of the bulk cannot be ad- 
dressed directly. Such questions include for example the fluctuations of the bulk density around p c , 
the fluctuations of the maximum around (p — p c )L, or the size of the second largest component. 

The contribution of this paper is a strong version of the equivalence of supercritical ensembles 
that provides a complete description for the thermodynamic limit and justifies the aforementioned 
picture. Precisely, if r\ G Xa is a configuration of particles on A we define 

M L (rj) = maxr] x 

and let m^fy) = argmax(7y) be the position where the maximum occurs. We can always enumerate 
the sites of A = {x±, . . . , xl} and define m^rf) to be the site with the smallest index should the 
maximum occur more than once. 
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We also define 

!Vx if x ^ y, z, 
r] y if x = z, 
i] z if x = y, 

and the operator T : Xa — ► Xa with Try = o- XL ' mL ^r] that exchanges the last and the maximum 
component of rj. 

We are ready to state the main result. 

Theorem 1. Let Tl be the a -field generated by r\ Xl , . . . rj XL . If p > p c , then 



lim^ sup \p N ' L o T- 1 [A] - v^- r [A] \ = 0. 

N/L^p 



This extends the result of Ferrari, Landim and Sisko [FLSj to the case where the number of sites 
increases to infinity together with the number of particles, and that of GroBkinsky, Schutz and 
Spohn |GSS| in the sense that convergence to the grand canonical distribution is obtained for the 
joint distribution under fi ' of all the components in the bulk. 

Given a measure p defined on a cr-algebra B, let || • || t . v . stand for the total variation norm 

\\p\\t.v. = sup \n(A)\. 

It is not hard to see that Theorem [T] then implies that 

^0, 



N,L _ ±_ \ "* ,,N,L „ „x,x L 

L 



xeA 

where v N,L is a probability measure on Xa with marginal on Tl-i given by v^ -1 , and such that 
the distribution of rj XL given !Fl-i equals the Dirac measure at N — ^2f = i Vxj- 

Several interesting facts about the invariant measures of the zero range process at supercritical 
densities are now simple consequences of Theorem [T] In view of ([2]) we would like to compute the 
fluctuations of Ml(ji) around (p — p c )L. This question was raised already in [JMPj and has been 
numerically investigated by Godreche and Luck (see appendix A. 2. 2 in |GL| ). The numerical ex- 
periments suggest that for b > 3 the fluctuations of Ml are of order yfh and Gaussian, while for 
2 < b < 3 they are of order . Theorem [T] and the obvious equality 

L-l 

M L (r,)=N-^2(Tr ) ) x , fi N ' L -a.s, 

x=l 

imply that the fluctuations of the maximum component reduce to the fluctuations of the sum of 
L— 1 independent random variables with mean p c around p c (L—l), for which standard central limit 
theorems are available 1GK] , The precise result is the following: 

Corollary 1. Suppose p > p c . 

a) If b > 3, that is if has finite variance a 2 = ^l^y^b-i) > ^hen for all x G M: 



lim u 

N/L^p 



N.L 



M L (rj) - (N - p c L) 
aUT2 <* 



/27T 



2 ' 2 du. 
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b) If b = 3, then for all x £ 



lim u 

N/L^p 



N,L 



M L {rj) - (N - p c L) 
2ViIogi 



< x 



2 / 2 du. 



c)If2<b<3, then for all x € 



lim u 

N/L^p 



N.L 



M L (n) - (N - p c L) 
(T(b)L)^ 



< x 



where C a is the density of the completely asymmetric stable law with index a and characteristic 
function tp(t) given by: 



logV(t) 



f - 1 - it*) J^T = ~C a \t\ a (l + i sgn(t) tan ^) 



Note that for b = 3 we still have Gaussian fluctuations after proper scaling. 

Clearly, one can go on and obtain limit theorems for the statistics of any order under p N,L from the 
corresponding result for product measures. For instance, the second largest component is given by 

M?\r,)= max (Tr,) x 

and the following limit theorem is a direct consequence of Theorem [T] and the estimate ([7]) for the 
tail probabilities under v^ c . 

Corollary 2. Suppose b > 2 and let p > p c . Then, for any x > 



lim a 

N/ L — > p 



N.L 



M L 2 \ V )< x(T(b)L) b -' 



The fluctuations of the bulk are closely related to the fluctuations of the maximum component. It 
follows from Corollary [5] that in the limit, is the only site where the number of particles is of 
order L. Given £ g (0, p — p c ) we define the bulk configuration as r/* = ^aT^^L}, an d the rescaled 
bulk fluctuation process Yl(-) € -D[0, 1] as 

, [Lt] 

YL(t) = —J2(v* Xj -Pc), 



where 



[ o-Vl 



ol 



{ 2/LfogT 
[(T(b)L)^ 



if b > 3 
if 6 = 3 
if 2 < b < 3. 



(3) 



The following corollary follows easily from Theorem [T] and Donsker's invariance principle or its 
extension by Skorokhod (Theorem 2.7 in |Sk| ) to i.i.d. random variables in the domain of attraction 
of a stable law. 
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Corollary 3. Suppose p > p c and let b > 2. Then under p N ' L 

Y L (-) &(•), as N -> oo, L — > oo, N/L -> p, 

where ^ is a standard Wiener process ij b > 3, or a completely asymmetric stable process with index 
a = b — 1 and characteristic exponent 

logVK-t) = jT - 1 - Hx) = -C a \t\ a (l - i sgn(t) tan ™), 

i/2 < 6 < 3. 

It is worth comparing Corollary [3] with the bulk fluctuations at criticality. If N — [p c L] then, 
according to a result in Thomas Liggett's dissertation (cf. Theorem 4 in [L]), Yl(-) converges in 
distribution to the bridge of conditioned to return to the origin at time 1. 

Theorem Q] can be also applied to the numerical simulation of the invariant states p N ' L , when 
N/L — ► p > p c . For large L, instead of drawing a sample from a distribution p N ' L , it is computa- 
tionally more efficient to draw L — 1 independent samples from a distribution v§ c , and assign the 
rest of the mass to a site uniformly distributed in {1, 2, . . . , L}. 

We present the proof to the main result in the following section. We conclude (section 4) by 
discussing two questions that arise naturally from Theorem [T] In the first one, we study a model 
such that the associated invariant measure v^ e has a stretched exponential tail, and prove that 
Theorem [T] still holds. In the second one, we consider a family of systems with particle numbers N 
deviating moderately from the typical value p c L, and refine our estimate of the threshold of values 
for N where a phase transition occurs. 



3 Proof of Theorem [T] 

We begin this section with a few observations on the model. Recall from the previous section that 
the jump rates are given by g(k) = 1 + r for k > 0, and the critical fugacity <j> c is equal to 1. Recall 
also that since we assume b > 2 both Z(cj> c ) and p c are finite. Although the precise value of the 
partition function, the critical density, or other statistics of are not important, it was pointed 
out in [GSS] that they can be explicitly computed using the hypergeometric identity pi] 

^ T(u + k)T(v + k) _ T(u)T(v)T(w -u-v) 

^ T(w + k) kl T(w - u)T{w - v) ' ^ 

valid for any u, v, w > with w > u + v. For instance, 

,t J i -P 1-. o j-l— f-.-^ —2 ( ) 



-, and if b > 3 then a 



6-1' r 6-2' (fe-2) 2 (6-3) 
We will next derive a smoothness estimate for the function 

W(k)-^ In -fcl- 1 _ (b-l)T(b)kl 

w{k) - „, e [Vx - k \ - z{ ^ )g{k)[ - r(fc + 6 + 1) ■ 

It is clear that W is decreasing, while from the elementary inequality 

1 + x > e^fe x > — 1 , (5) 
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one can easily deduce that W(k)k b is increasing. Thus, for k\ < k2 we get 



, \ b 

«i 



W(ki) > W(k 2 ) > W{kx) (6) 



We can also apply (U|) to compute the tail probabilities of v^ c as follows 



£ W(k) = (b-l)T(b)Y^ 



r(m+l + fc) _ T(b) ml 



r(m + b+l + k) Ttm + b) 

k=m fc— / \ / 

Hence, if we denote by F the distribution function of and by F — 1 — F its tail, we get the 
following asymptotic behavior at infinity 

W(k) ~ (6 - l)r(6)fc- fc , and F(x) ~ T{b)x 1 - b . (7) 

This observation explains the normalizing constants in the statements of Corollary [2] and Corollaries 
[T] and [3] for 2 < b < 3. The logarithmic correction when 6 = 3 comes from the direct computation 

L 1 

E"* H Mv*<L}] ~4^-~41ogL. 

fe=i 

The proof of Theorem [1] relies on a local limit theorem for the (unconditioned) measure at criticality. 
It estimates the probability of the event we are conditioning upon in the definition of p N,L . Such 
a result first appeared in Nagaev |N1] for b > 3 and Tkacuk [T] for b < 3. Baltrunas [B] gives an 
accessible proof that encompasses all values of b > 2. 



Proposition 1. If p > p c then 



N,L 



Hm ■jfrM-fl ( 8 ) 



Equation ([5]) says that the most probable way that the rare event {Sl(t]) — N} occurs is when one 
variable takes up all the "excess mass", while the remaining L — 1 ones assume typical values. This 
behaviour is to be contrasted with the large deviations behavior for random variables with finite 
exponential moments, where the rare event is realised by all variables taking values close to the 
atypical p. 

We proceed now with the proof of Theorem [T] 

Proof of Theorem]^ Recall from Section 2 that a Xi,Xj stands for the mapping that exchanges the 
i—th and the j-th components of 77, and that T denotes the transformation that exchanges the last 
and the maximum components of 77. 

Let A C {7/ : r] XL > i] Xj , j — 1,2, . . . ,L — 1}. Due to the invariance of p N,L under a XL,Xi , I = 
1, . . . , L, we get 

L 

» N < L [T- X A] = ^ N ' L [ T ~ lA n {m L = £}] 

e=i 

L 

t=\ 

_ v^ c [A n {Sl(v) = N}] 

ul[S L ( V ) = N] ■ W 
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Consider a sequence Cl such that Cl/L — > and Cl/o-l — > oo, where is defined in ([3]). Let 
D L = {m: \N — p c L — m\ < C L }, t L = N - p c L - C L and B L = {r/ : r] XL e D L ; max\<j<L-i Vxj < 
t L }- 

Suppose now that A E Tl-\ = &{t) x1 , . . . ,rj XL _ 1 } We will apply equation (J9j) to the set A fl Bl- 
Note that 

i£ [A n Bl n {5 L (?7) = TV}] = 51 Wim)^- 1 \a.C\ { ^ 7j x — N - m; max?], < t L } . 

mED L x 

In view of ((6|) we can replace each value W^m) in the range of summation by W(N — [p c L]), creating 
an error that is negligible as L — > oo uniformly in A. That is, 



W(N - [ Pc L]) v\~ x A n || Vx - PcL\ < C L - maxrfa <t L } 



(l + o(l)). 



Since Cl/o-l — * oo, the central limit theorem implies that 



.L— 



'[IE 



Vx - PcL\ < C L 



1 as L — > oo, 



and there is also the elementary estimate 

i/^^rnaxryx < t £ ] = (l - F{t L )Y * — >1 as L -> oo. 
Combining these two observations, we get 

4 c [A n Bl n = JV}] = I^(iV - [p c i]) z^ 1 [A] (l + o(l)) . 

Together with equation ([9]) and Proposition [1] this establishes that 



lim sup 



N/i — p 



^ L oT-i[AnB L ]-v^[A] 



= 0. 



In particular, if A = Xa, we get that 



N/L~,p 

B C L = Xa \ Bl- The assertion of the Theorem now follows by combining the last two equations. □ 



4 Remarks 

We identified the condensation phenomenon present in supercritical zero range processes by proving 
the equivalence of ensembles in the standard Evans' model. It should be clear however that the 
essential ingredient for the proof is a Local Limit Theorem in the form of {5} . There are thus two 
possible directions to generalise Theorem [T] Its validity should be established for a greater variety 
of models, and the point where the phase transition with the emergence of a large cluster occurs 
should be determined with greater accuracy. 

We describe next how the proof can be adapted to a model for condensation with stretched ex- 
ponential tails, also proposed by Evans. 



9 



Suppose the jump rates are given by the function g with 

m = { l + l iik>0 do) 

JV [0 if fc = 0, 

where A £ (|, 1). The critical fugacity is still 1, although it is not possible to explicitly compute the 
distribution function and the critical density in this case. Nevertheless, it is elementary to see that 
W(k) = \r) x = k\ is decreasing while W(k) exp ) is increasing so that we have 

W(k 1 )>W(k 2 )>W(k 1 )exp^-/3 2 _ 1 j, fe x < fc 2 . (11) 
In fact, using © one can check that 



*(Jb)<exp(-£^)<Ce^(-£ 

\ m=l / 



-A 



and the following asymptotic behavior for W holds 

W(k) <~ A exp ^ — —J as k — > oo. 

This yields the asymptotic behavior of F(x) 

W ~ — — exp I — as x — > oo. (12) 

p V 1 — A / 

In this context, Nagaev |N2| has proved that © is satisfied as long as N = p c L + ^f{L)L^ with 
j(L) — > oo as L — > oo. In view of equations (fTT|) and (fl2|) . we may choose the sequence Cl = V L log L 
in the line following the expression ([5]), and adapt the arguments presented in the previous section 
to prove the following theorem. 

Theorem la If g(-) is given by l[T(J\) and N = p c L + j(L)L"T\ where lim7(L) = oo, then 

lim sup \p N ' L o T- 1 \A] - v L r x \A] I = 0. 

In a similar fashion we can relax the conditions on N in Theorem [1] provided we prove the validity 
of ijSJl for values of N deviating only moderately from its typical value. For instance, when b > 3 
Theorem 2 in [Do implies that if (N — p c L) j\fL — ► oo then 

i£ [S L (v) = N]=-j=<p (^p) (1 + o(l)) + LW(N - [ Pc L]) (1 + o(l)), 

where (p(-) is the density of the standard normal distribution. It is not hard to see that in this case 
(IH) holds as long as 

N = p c L + ^v^gL (l + + , with li^(L) = oo. (13) 

Once again, choosing Cl = \/L logi we can prove the following refinement of Theorem [T] 

Theorem lb If g(-) is given by (Qp with b > 3 and N is as in U3\) , then 

lim sup \p N - L o T- 1 [A] - v L r x [A] I = 0. 
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Similar refinements of Theorem [T] can be obtained for the case when b < 3. 
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